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A novel method of solving scattering problems for bound pairs on a lattice is developed. Two 
different break ups of the hamiltonian are employed to calculate the full Green operator and the 
wave function of the scattered pair. The calculation converges exponentially in the number of basis 
states used to represent the non-translation invariant part of the Green operator. The method is 
general and applicable to a variety of scattering and tunneling problems. As the first application, 
the problem of pair tunneling through a weak link on a one-dimensional lattice is solved. It is found 
that at momenta close to ±7r the pair tunnels much easier than one particle, with the transmission 
coefficient approaching unity. This anomalously high transmission is a consequence of the existence 
of a two-body resonant state localized at the weak link. 
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Introduction. Scattering of bound particle complexes 
has been a major subject of atomic, molecular and nu- 
clear physics for decades. In "lattice" solid state physics 
the prime system of interest has been the excitoni^, in 
which the constituent particles, an electron and a hole, 
have different masses. The bound pair of two magnons 
in lattice magnetism is an example of a complex with 
equal masses 3 . In recent years, the concepts of lattice 
bipolarona^*Si& and bisolitons^ have been developed, 
in particular in relation with high-temperature supercon- 
ductors and conducting polymers. 

Many properties of these particles derive just from 
their composite nature rather than from the particu- 
lars of the binding interaction. They can therefore be 
studied within the framework of the "generic" two-body 
system, in which a model potential is introduced to en- 
sure binding, yet the simplicity of the potential enables 
rigorous analysis of the quantum mechanical problem. 
This approach has been popular and the physics of two- 
particle bound complexes in translation invariant lattices 
is now well-understood, see for example Refs. [T(illlllll2| 
and the bibliography therein. Much less is known about 
non-translation invariant cases. When defects or bound- 
aries are present the two-body problem can no longer 
be reduced to a one-body problem, which significantly 
complicates analysis. In continuum physics, scattering 
of bound pairs was approached from the general three- 
body formalismiii^ii^, although no exact results were ob- 
tained beyond the one dimension with delta-function po- 
tentials. On a lattice, the previous research was limited 
to the surface excitonic effectsi£*i£. Bulato\*i£*i& devel- 
oped a general theory and an efficient numerical proce- 
dure to obtain the energy spectra and wave functions of 
lattice excitons in the presence of a surface. 

In this Letter, we extend the method of Refs. 
to the general scattering problem of lattice bound pairs. 
The method consists of calculating the full two-particle 
Green operator G and then acting with it on the wave 
function of an incident pair ^y. The core feature of 
the method is the usage of two different decompositions 



of the hamiltonian on a zero part and a perturbation. 
The first decomposition is applied to find G while the 
second decomposition is used to calculate the scattering 
amplitudes. The accuracy of the method increases ex- 
ponentially with the number of lattice sites used to ap- 
proximate the non-translation invariant part of G. As 
the first application of the method we solve the problem 
of tunneling of a one-dimensional bound pair through a 
weak link on a chain. We find that the pair transmis- 
sion at large lattice momenta is significantly enhanced in 
comparison with the transmission of a single particle. In 
fact, the transmission coefficient approaches unity at the 
Brillouin zone boundary. 

Method. The generic model consists of free motion of 
two particles Hq, interparticle interaction V (which is 
usually attractive), and single-particle scattering U: 

H = H a + V + U (1) 
= H v + U (2) 
= Hu + V. (3) 

Here two partial hamiltonians Hy = Hq + V and Hjj = 
Hq + U are introduced. Equations J5J) and (0 define 
the two decompositions mentioned above. Using the de- 
composition the- full wave function \I/ satisfies the 
Schrodinger equation 

* = *y + Gy[/*, (4) 

where Hy^y — E^y, &v has the appropriate boundary 
conditions at infinity, and Gy{E) = (E — Hy + ij)^ 1 is 
the Green operator of Hy. Three other Green operators 
G, Go, and G\j are defined analogously. In the basis of 
localized lattice states the Green operators can be repre- 
sented as ordinary matrices, albeit of infinite size. 

Ordinarily, equations like J3J are used to develop per- 
turbative expansions for from the knowledge of the 
partial Green operator Gy. Now suppose that the full 
Green operator G is known. Since G = (1 — GyU)~ 1 Gy 
and * = (1 - G v U)~ ll $>v, the last term in Q is re- 
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arranged as follows 

G V U^ = GvUil-GyUy^v 

= (1 - GvUy-GvUVv = GU^v, (5) 

so that, from Eq. Q 

* = tfy + GU^! V = (1 + GUpv (6) 

Thus if G is known, the full wave function can be found 
from the last equation by matrix multiplication. 

Now comes an important observation. Since G is the 
full Green operator it does not matter how it is obtained. 
In particular, one is not obligated to use the same de- 
composition (J2J) that has led to Eq. ©. For scattering of 
bound pairs it is more convenient to use the decomposi- 
tion Q, which yields 

G=(1-G U V)- 1 G U = A- 1 G U . (7) 

The advantage of this approach is that Gjj is the Green 
operator of two non-interacting particles, both scattered 
off the potential U. Therefore Gu can be calculated as a 
convolution of two one-particle Green operators gjj: 

f°° de 

G u (r 1 r' 1 ;r 2 r' 2 ;E) =i —g u {r 1 r 1 ;e)g u (r- 2 r' 2 ;E - e). 

(8) 

In turn, gu follows from solving of a one-particle scatter- 
ing problem: 

gu = (1 - .got/r^o, (9) 

where go = (E — Hq + 17) is the one-particle Green 
operator of the translation-invariant system. The zero 
operator g is most easily calculated from the spectral 
expansion 

gik(r-r') 

where £k is the one-particle spectrum. Gjj can also be 
calculated from the two-particle spectral expansion^. 
Thus the strategy of the present method is to use the 
decomposition © and formulas l(7|)- (|10|) to obtain the 
Green operator G, and then use the decomposition 
and formula JJJJ to calculate the full wave function \& and 
the scattering coefficients of interest. 

Calculation of (1 — GjjV)^ 1 . Once Gjj is known, the 
main task is to invert the matrix A = 1 — GjjV, see 
Eq. J7J. The way of calculating A^ 1 is the second key 
component of the present method. Observe that invert- 
ing (1 — GjjV) is analogous to inverting (E — H), i.e., to 
calculating a Green operator. Imagine a Gjj that con- 
sists of a translation-invariant part G\j and a perturba- 
tion SG = Gij — Gy which is localized in real space. Then 
the translation invariant part (1 — GyV) plays the role 
of the translation invariant part of (E — H) while 5GV 



the role of the localized perturbation. Performing the 
standard transformation one obtains 

A' 1 = (1 - G%V - SG V)- 1 = (1 - B^SG V)~ l S -1 , 

(11) 

B = 1 - G V V. (12) 

Thus inversion of A is replaced with two inversions. 
The first inversion is that of B. Since B involves only 
translation-invariant matrices this is achieved by chang- 
ing to the quasi-momentum representation in which B 
is block-diagonal with the block size equal to the range 
of V— in relative coordinates. The second inversion is 
that of (1 - B^SG V). The latter is the sum of the 
unit matrix and a matrix localized around the scatter- 
ing region, which is due to the localization of SG. Thus 
only inversion of a finite-size matrix that contains the 
non-zero elements of B~ 1 8G V is required. As a result, 
the abstract task of inverting the infinite matrix A is re- 
placed with two easy-to-perform operations on finite size 
matrices. 

To summarize, the algorithm begins with the calcula- 
tion of the Green matrix Gu from Eqs. jSjl- lfTUI) . Then 
Gjj is separated into the translation-invariant part G v 
and the remainder 5G. On the next step, A^ 1 is calcu- 
lated according to Eq. (|llfl . and then the full Green op- 
erator is obtained from Eq. Q. Finally, an unperturbed 
pair wave function \Ev is chosen and the scattered wave 
function is calculated from Eq. This formulation is 
completely general, and can be applied to a variety of 
particular cases. One such problem is analyzed below. 

A chain with a weak link. Consider a one-dimensional 
chain characterized by the nearest-neighbor hopping ma- 
trix element t > and the Hubbard-like attraction of 
strength \v\ > 0. The hopping amplitude between sites 
n = and n = 1 contains an additional element t' . The 
resulting hamiltonian reads 

H = -t ^ C laCn'a ~ M^ C «T Cn T C il C «i 
(nn')cr n 

+t'^(4 CT Cla + cl £r C 0CT ), (13) 
a 

where (nn'} denotes pairs of nearest neighbor sites. The 
value t 1 = corresponds to the absence of any scattering 
while t' = t corresponds to two decoupled semi-infinite 
chains. A standard solution of the one-particle scattering 
problems yields the transmission coefficient: 

(t'/t-l)(e ik -e- ik ) 
k e~ ik - (t'/t - l) 2 e lk ' V ; 

where k is the one-particle momentum. The modulus of 
the transmission coefficient \t^ \ is shown in Fig.^ m bold 
lines. Note that = at k = or k = tt. 

In the absence of scattering (f' = 0), two particles 
form a singlet bound state with an (unnormalized) wave 
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FIG. 1: Bold lines: modulus of the single particle transmis- 
sion coefficient \rk \ through the weak link, see Eq. 11411 . From 
the top curve down: t' /t — 0.1, 0.2, . . . , 0.9. Thin lines: the 
same quantity in the presence of a resonant state at the top 
of the single-particle band, see Eq. 1181 . 



function 



^ v {n u n 2 ) = e 



_ ±i^-(n 1 +n 2 ) -\\n!-n 2 \ 



(15) 



where K > is the total momentum of the pair and 
sinhA = |«|/[4t cos (K/2)]. The energy of the bound 
state is E — -y/v 2 + 16t 2 cos 2 (K/2) < 0. We choose to 
study scattering of pairs incident from the left with en- 
ergy E < —At to prevent the processes of pair breaking in 
two free particles. At these energies the full wave function 
© has the asymptotic '5 — > ^y + R^f v at n l7 n 2 — * — oo, 
and 'F — ► T^ v at n,x,n 2 — > +oo. We are interested in 
the pair transmission coefficient T as a function of the 
pair momentum K and model parameters t' and \v\. 

Determination of T begins with calculating Gjj from 
Eqs. JSJl-JTSJ using as input et = —2tcosk and U that 
has all the matrix elements zero except uoi — Uio = 
t'. The translation-invariant part of Gjj can be obtained 
numerically by simply setting t' = 0. Alternatively, the 
two-particle spectral expansion yields for at E < —At 
the following expression [for the model (|13fl . only n\ = 



n 2 , n x 



= nl matrix elements of G? r are needed because 



of the locality of the Hubbard attraction] : 



G;j( n, n ;n\n / ) 



dkidk 2 cos ki(n — n') cos k 2 (n — n') 



„ (2n) 2 E + 2t cosh +2t cos k 2 
T dq cos q(n — n') 



2tt ^Je 2 - 8t 2 - 8t 2 cosq 



(16) 



In accordance with the general scheme, the matrix 6G is 
calculated by subtracting G^ from G\j. SG(n, n') is very 
localized around the weak link n,n' = 0, 1, see Fig. [21 

The last thing we need is an expression for B^ 1 , see 
Eq. (|12fl . Again, only the matrix elements in the block 




FIG. 2: The non-translation invariant part of the Green op- 
erator 5G(n,n) for \v\ = 4.1 1 and t = 0.3 1. The pair mo- 
mentum is K — 0.3-7T (left panel) and K = 0.8-7T (right panel). 
The radius of the circle represents the modulus \SG\. Notice 
the high degree of localization around the weak link. 



1 



incident and reflected wave 



Wanstnitted wave 



FIG. 3: The pair wave function ^ k (tii, n-2) for K = 0.3 7r, 
\v\ = 4.1 1 and t' = 0.3 1. The radius of the circle represents 
the modulus |\E'|. 



711 = n 2 and n! x = n' 2 are required. By diagonalizing the 
block by a Fourier transformation one can show thati&Sa 



B- l (n;n') 



* dq 



cos q(n — n') 



— 7t 2?r 1 



(17) 



^/(£+ 47 ) 2 -16t 2 cos 2 (g/2) 



Numerical results. The results obtained in the pre- 
ceding section enable calculation of the full two-particle 
Green operator, the exact pair wave function, and the 
scattering coefficients of bound pairs for the model i|13[l . 
In Fig. |3|we show the modulus \^ k\ as a function of the 
particle coordinates n± and n 2 for K = 0.3 7r, \v\ = 4.1 1 
and t' = 0.3 1. Notice how reflection off the weak link cre- 
ates interference between the incident and reflected wave 
functions. In contrast, the transmitted wave (in the lower 
right part of the graph) has a constant amplitude. 

In Fig. 01 we show the pair transmission coefficient Tr-. 
As a function of pair momentum, Tk behaves qualita- 
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FIG. 4: Pair transmission coefficient Tk for \v\ = 4.1 and dif- 
ferent t' It. From the top curve down: t'/t = 0.1, 0.2, . . . , 0.9. 

tively different from the one-particle transmission Tk , see 
Fig. H r fc first increases with momentum but then de- 
creases and vanishes at fc = n. In contrast, Tk is a mono- 
tonically increasing function of momentum, and reaches 
unity at K = %. Thus at large lattice momenta a bound 
pair is transmitted through a weak link much easier than a 
single particle. The likely physical reason for the anoma- 
lously high transmission is resonant tunneling through a 
two-body state localized at the weak link. Bulatov and 
Danilo\i2i previously analyzed the two-particle spectrum 
of a semi-infinite Hubbard chain, i.e. model 1)13(1 with 
t' = t. They found that the chain boundary introduces 
a resonant state with E = — |i>|, i.e. exactly at the top 
edge of the pair band. We conjecture that such a state 
exists also at t' ^ t and facilitates efficient transmission 
through the weak link of pairs with energies close to the 
top of the band, i.e. with momenta close to tt. 

It is instructive to compare this effect with one-particle 
tunneling through the weak link in the presence of a res- 
onant state. Such a state appears in the model (| 1 31) with 



\v\ = if an additional one-particle repulsive potential 
w is added at the two sites on either side of the weak 
link. At id = t', the state has the energy of the top of 
the one-particle band, E — 2t. For those parameters, the 
transmission coefficient 



(e ik ~e- lk )-2(t'/t)(l + e lk )' v ; 

This function is shown in Fig^i 11 thin lines. The reso- 
nant state qualitatively changes the transmission at large 
momenta. Instead of vanishing ?k actually approaches 
unity. The overall shape of the curves is remarkably sim- 
ilar to pair transmission curves of Fig. 01 which further 
supports our interpretation of pair tunneling as through 
a resonant state. 

Summary. We have developed an efficient procedure of 
calculating scattering coefficients of bound pairs on a lat- 
tice. The key technical advance of the paper is the usage 
of two different decompositions of the hamiltonian; one 
is used to calculate the full Green operator of the system 
while another to find the resulting wave function of the 
pair. Another important element is the method of invert- 
ing the matrix (1 — GjjV), which is based on separating 
Gjj on a translation invariant part and a part localized 
around the scatterer, see Eqs. ijTTJ and (|12fl . The numer- 
ical accuracy of the method scales exponentially in the 
number of basis states chosen to represent the localized 
part. As formulated, the method is quite general enabling 
accurate investigation of a variety of scattering and tun- 
neling problems. As the first application, we have studied 
transmission of bound pairs through a weak link on the 
one-dimensional chain. Contrary to simplistic expecta- 
tions, we have found that at large momenta the pairs 
penetrate the barrier easier than single particles. The 
anomalously high transmission has been identified with 
tunneling through a resonant pair state. More two-body 
scattering problems are currently under investigation. 
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